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ABSTRACT 

A series of high-resolution numerical simulations were performed to study the structure and sub- 
structure of Milky Way (MW)- and cluster-sized halos in a A— Cold Dark Matter (CDM) cosmology 
with self-interacting (SI) dark particles. The cross section per unit of particle mass has the form 
o~dm = 0o(1/vioo) Q j where (To is a constant in units of cm 2 gr~ 1 and vioo is the relative velocity in 
units of 100 kms -1 . Different values for <tq with a — or 1 were used. For small values of cfdm= 
const. (<; 0.5, a = 0), the core density of the halos at z = is typically higher at a given mass for 
lower values of Co or, at a given a , for lower masses. For values of ao as high as 3.0, both cluster- and 
MW-sized halos may undergo the gravothermal catastrophe before z = 0. The core expansion occurs in a 
stable regime because the heat capacity, C , is positive in the center. After the maximum expansion, the 
isothermal core is hotter than the periphery and C < 0. Then, the gravothermal catastrophe triggers. 
The instability onset can be delayed by both the dynamical heating of the halo by major mergers and 
the interaction of cool particles with the hot environment of a host halo. When a — 1, the core density 
of cluster- and MW-sized halos is similar. Using cfdm = 0.5 — 1.0 (1/vioo), our predictions agree with 
the central densities and the core scaling laws of halos both inferred from the observations of dwarf 
and low surface brightness galaxies and clusters of galaxies. Regarding the cumulative v max — function of 
subhalos within MW-sized halos, when (o- n ,a) = (0.1,0.0), (0.5,0.0) or (0.5,1.0) it agrees roughly with 
observations (luminous satellites) for v max > SOkms^ 1 , while at v max = 20kms _1 the functions are al- 
ready a factor 5-8 higher, similar to the CDM predictions. For (cto, a) = (1.0,1.0), this function lies above 
the corresponding CDM function. The structure and number of subhalos are affected by the scattering 
properties of the host halo rather than by those of the subhalos. The halos with SI have more specific 
angular momentum at a given mass shell and are rounder than their CDM counterparts. However, the 
angular momentum excess w.r.t. CDM is small. We conclude that the introduction of SI particles with 
o~dm oc 1/vioo rnay remedy the cuspy core problem of the CDM cosmogony, at the same time keeping a 
subhalo population similar to that of the CDM halos. 



1. INTRODUCTION 

The current paradigm for cosmic structure formation 
is based on the cold dark matter (CDM) cosmological 
models, where CDM particles are assumed to be collision- 
less. This is the simplest assumption about the nature 
of these particles, not yet detected. The predictions of 
the most popular encarnation of the CDM variants, the 
spatially flat with a non-vanishing cosmological constant 
model, are in excellent agreement with a large body of 
observational data at large scales. With the advent of 
high-resolution numerical simulations and new observa- 
tional techniques, a comparison between models and ob- 
servations at small scales has become possible. Apparent 
conflicts have emerged as a result of this comparison: (i) 
the predicted halo inner density profiles are cuspy, in dis- 
agreement with the shallow cores favored by observations 
of dwarf and low surface brightness galaxies (LSB), and (ii) 



the predicted number of subhalos within Milky Way-sized 
(MW-sized) halos overwhelm the observed abundance of 
satellite galaxies in the Local Group. 

Recently, a plethora of alternative theories, which mod- 
ifiy the predictions of the CDM model at small scales but 
retain its successes at large scales, have been proposed 
(e.g., Dave et al. 2001 and references there included). In 
one of these scenarios, CDM particles are assumed to be 
self-interacting (SI) in such a way that the heat flux to 
the core smooth out the density cusp, and simultaneously 
reduces the amount of substructure by evaporating orbit- 
ing subhalos (Spergel & Steinhardt 2000). Several authors 
have explored this model numerically and analytically and 
concluded that the relevant regime for structure forma- 
tion would have to be the optically thin one (Moore et 
al. 2000; Yoshida et al. 2000a; Firmani et al. 2001b; 
Kochanek & White 2000; Dave et al. 2001; Hennawi & 
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Ostriker 2001; Balbcrg, Shapiro, & Inagaki 2002). By 
means of cosmological N-body simulations, using a con- 
stant cross section per unit mass, a dm, Dave et al. (2001) 
found that halos at galaxy scales have long-lived shallow 
cores whose sizes agree with observational inferences when 
o~dm ~ 5 cm 2 gr _1 . At the cluster scales, Yoshida et al. 
(2000b) also found that the SIDM halos present long-lived 
shallow cores. According to them, odm should be smaller 
than w 5 cm 2 gr _1 in order for their results to be in agree- 
ment with observations of cluster of galaxies. 

On the other hand, N-body simulations of relaxed, iso- 
lated, SI dark matter (SIDM) halos with initial cuspy den- 
sity profiles and with a relatively high a dm value show a 
core evolution which is remarkably rapid. After the maxi- 
mum soft core radius is reached, in time scales of the order 
of the halo dynamical time, the core collapse is triggered 
almost immediately (Kochanek & White 2000, hereafter 
KW; Burkcrt 2000; see also Quinlan 1996). Dave et al. 
(2001) argue that the fast core collapse seen in these sim- 
ulations may be attributed to the absence, in the isolated 
case, of infalling dynamically-hot material. Dave et al. 
also find that the central density of simulated cosmologi- 
cal halos (at 1 h _1 kpc from the center) remains approxi- 
mately constant with halo mass, a result actually not ex- 
pected for a constant a dm (Yoshida et al. 2000b; Firmani 
et al. 2001a,b). In fact, Firmani et al. (2001b) find that 
p c is nearly independent of halo mass only when a dm is 
inversely proportional to the relative velocity of colliding 
particles. They used a numerical code based on the colli- 
sional Boltzman equation and took into account the cos- 
mological mass aggregation process during the formation 
of their spherical symmetric halos. They obtain a central 
density of p c w 0.02M Q pc~ 3 (a value inferred from ob- 
servations, h = 0.65 was assumed, Firmani et al. 2000a) 
for odm ~ (0.6/vioo) cm 2 gr _1 , where vi o is the relative 
velocity vi 2 in units of 100 kms -1 . 

As suggested by all these studies, important questions 
regarding the SIDM cosmology still remain without a sat- 
isfactory answer or without an answer at all. Do the SIDM 
cosmological halos undergo a core collapse on scales times 
shorter than a Hubble time? How does the cosmological 
mass aggregation history (MAH) affect the core evolution? 
Does a constant cross section produce halos with a nearly 
constant central density regardless of their mass? Which 
are the predictions for SIDM models with a cross section 
inversely proportional to V12? How is the subhalo popu- 
lation at cluster and galaxy scales in a SIDM cosmology? 
In this paper, we will address these questions using high- 
resolution cosmological simulations of cluster- and galaxy- 
sized halos with a dm °c const and oc l/vi 2 . 

Our study is aimed to explore whether the SIDM model 
predictions are in better agreement with observations than 
the CDM ones. As already disscused above, the key pa- 
rameter for SIDM is the interaction cross section, 

VDM = O- (^) • (1) 

A number of authors have attempted to constrain the 
range of values for the pair of parameters (cro,a) (Fir- 
mani et al. 2000, 2001a; Miralda-Escude 2000; Wandelt 
et al. 2000; Meneghetti et al. 2001; Gnedin & Ostriker 

1 In ART, only those cells that are not refined are allowed to own 



2001; Hennawi & Ostriker 2001). According to analyti- 
cal and semi-analytical arguments of Hennawi & Ostriker 
(2001; HO hereafter), a dm should be large enough to pro- 
duce a shallow core in agreement with observations, but 
small enough so the core halo would not collapse in a Hub- 
ble time or subhalo evaporation within cluster-sized halos 
would not alter dramatically the elliptical Fundamental 
Plane relations. This leaves the SIDM cosmology with a 
limited range of ((To, a) values. HO also introduce an ex- 
tra strong constrain for SIDM, the mass of a central black 
hole in the MW and the lack of it in M33. According to 
them, accretion of SIDM onto seed black holes produces 
holes that can be much more massive than those observed 
if a dm is roughly as large as needed to form the soft cores 
inferred from observations. In §5.2 we will discuss this 
constrain in light of our results. 

In §2 we describe the method and strategy to simulate 
SIDM halos, and present the simulation of a monolith- 
ical halo (§2.1) as well as the cosmological simulations 
carried out in this paper (§2.2). In §3.1 we analyze the 
density profiles and central densities measured in our sim- 
ulated cluster- and galaxy-sized halos with different values 
of (<To, a) and different merging histories. In §3.2 we dis- 
cuss the concentration parameters of halos and subhalos, 
while in §3.3 and 3.4 we present results regarding the el- 
lipticity and angular momentum distribution of the SIDM 
halos. Section 4 is devoted to the study of the subhalo 
population in the MW- and cluster-sized SIDM halos. In 
§5 we discuss our results and the viability of the SIDM 
models, and in §6 we summarize this paper and present 
our conclusions. 

2. SIMULATIONS 

The Adaptive Refinement Tree (ART) N-body code 
(Kravtsov, Klypin, & Khokhlov 1997) has been used to 
run the N-body simulations. The ART code achieves high 
spatial resolution by refining the base uniform grid in all 
high-density regions with an automated refinement algo- 
rithm. The clastic collisions, characterized by a scattering 
cross section a dm, are implemented in ART in the follow- 
ing manner. In a medium with density p, the predicted 
distance d that a particle travels without collisioning is 
given by (Gibbs 1994) 

d=-AIn(l-P), (2) 

where A = 1/ pa dm is the mean free path and P is a 
random number distributed uniformly between and 1. 
A pair of particles collide if this distance d results lower 
than the distance vi^Ai, where V12 is the relative veloc- 
ity between the particle and one of its nearest neighbors, 
and At is the time step. In a inhomogenous medium, 
we substitute p by the local density at the particle's po- 
sition, pi. We exploit the mesh hierarchy structure of 
ART to evaluate both the local density and the relative 
velocity. On each level of the structure, starting from the 
finest level, ART assigns the density in each cell using the 
standard cloud-in-cell technique (Hockney & Eastwood 
1981). The density at the particle's position is found 
by interpolation 1 using cloud-in-cell. The search for a 
partner to collide with may involve many neighbor cells 
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and more than one level. We find that the search for 
partners of particles that lie in low-density regions fails 
often; this, however, does not affect our results because 
collisions in these regions are very rare. It happens many 
times that the partner is found in the same cell where 
the particle is, and thus the search is limited to only 
one cell. The new velocities after collision are computed 
by imposing energy and momentum conservation and 
by orienting the after-collision relative velocity randomly. 




Fig. 1. — Evolution of the density, 3D velocity disper- 
sion, and heat capacity profiles (top, central and bottom pan- 
els, respectively) for the NFW isolated SIDM halo of mass 

2 x 10 11 h -1 M0 and initial concentration cnpw = 14. Tha 
maximum expansion of the core occurs at about 2td yn . With 
a solid line is represented the analytical NFW function from 
which a relaxed and stable dark matter halo is reproduced 
(dotted line). We see that at 9 t dyn the central density has 
increased by almost a factor 2. 

2.1. Self- Interacting simulations of a monolithic halo 

In order to test our collisional algorithm on ART code 
we ran a simulation of a monolithic (isolated) halo of 

3 x 10 5 particles with an initial NFW density profile, and 
SI turned on. The virial mass of the halo is M vir = 
2 x 10 11 h _1 M Q , with a NFW concentration of c NF w = 14. 



If we assume the Hernquist model to have the same den- 
sity as the NFW one in the central r _1 region, we can 
relate the scale radius r# of the Hernquist density profile 
with cnfw- Given r# and a dm we can estimate the di- 
mensionless cross section a and the core relaxation time 
t rc , both as defined by Kochanek & White (2000, hereafter 
KW) who also simulated SIDM monolithic halos but with 
a Hernquist profile. The simulation was run when the SI 
algorithm was still on test and used the particle velocity 
for the computation of the collisional probability instead 
of V12; this produces a smaller scattering rate than when 
one uses v i2 (KW). The model has uum = 9.0 cm 2 gr _1 . 

The evolution of the density, 3D velocity dispersion, and 
heat capacity profiles are shown in Figure 1; time is mea- 
sured in units of the dynamical time, idyn, as defined in 
KW. The evolution proceeds in much the same way as in 
previous similar SI simulations (e.g., Burkcrt 2000; KW). 
If the initial equilibrium configuration of a halo has a NFW 
density profile, the velocity dispersion (temperature) pro- 
file peaks out of the center, close to the scale radius rNFW 
(dotted line in Fig. 1); i.e., there is an i nversion of the 
temperature profile. We have measured the potential, ki- 
netic, and total energy, W, T, and U, respectively, of our 
monolithic halo as a function of radius and found, for the 
initial configuration, that both the total energy (see also 
Fig. 2 of Lokas & Mamon 2001) and the heat capacity, C, 
of the innermost region are positive (see Fig. 1). 

Therefore, when the SI is turned on, heat starts flow- 
ing into the core from the immediate surroundings and 
the core expands. An increase in the core's total energy 
produces an increase in the core's temperature and in the 
ratio 2T/\W\; in fact, the maximum expansion of the core 
occurs approximately when 2T/\W\ also reaches its max- 
imum. At this point the core is isothermal and the tem- 
perature inversion is gone (Fig. 1; see also panel (a) in 
Fig. 2). Now, the isothermal core is hotter than the sur- 
roundings and heat starts flowing outwards (note that this 
process happens in the outer central region where C is al- 
ready negative). The central region, where both U and C 
are positive, reduces in size (Fig. 1). The core looses heat 
and contracts; the system is unstable and the gravothcrmal 
catastrophe phase triggers (see Binney & Tremaine 1987, 
§8.2): the core becomes hotter as it looses energy out- 
wards, and the increased temperature difference leads to 
an even faster energy loss with the consequent core shrink- 
ing (panel (b) in Fig. 2). 

An interesting question is whether the inverse instabil- 
ity, the gravothcrmal expansion, may arise sometime in 
the evolution of a halo. This bring us to the question: why 
at the begining of the SI evolution of the isolated halo, 
when there is an inversion of the temperature profile, a 
run-away core expansion does not occur? The reason is 
that the central region has positive total energy and heat 
capacity. Let us now analize the halo after the maximum 
core expansion and let us assume that for some reason (for 
example, kinetic energy injection by mergers), the outer 
region of a halo becomes hotter than the isothermal core 
(panel (c) in Fig. 2). The system in this case does not 
obey the conditions for the gravothermal expansion insta- 
bility because, again, U and C in the central region are 
positive (due to the strong halo heating, C may become 
positive even at the periphery). A flow of heat into the 
core makes it expand and brings most parts of the system 
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to an isothermal state (panel (c)), which might be unsta- 
ble to core collapse because the isothermal sphere becomes 
a local entropy minimum rather than a local entropy max- 
imum. Thus, the core collapse phase for cosmological 
SIDM halos can be delayed due to an external energy 
injection, but not reversed. We will return to this diss- 
cusion in §3 where we analyze cosmological SIDM halos. 



NFW (a) 




logR 



Fig. 2. — Schematical behavior of the 3D velocity disper- 
sion (temperature) profile in a SIDM halo, and of the loga- 
rithmic density profile (encapsulated boxes). Panels (a), (b), 
and (c) show the core expansion phase (stable), the gravother- 
mal catastrophe phase (unstable), and the core collapse delay 
(and core expansion) when the halo is externally heated, re- 
spectively. See the text for more explanations. C is the heat 
capacity. 



The notion of gravothermal catastrophe has been intro- 
duced for a non-singular isothermal model (Lynden-Bcll 
& Wood 1968). Although in the cosmological case the 
halos are not described by isothermal spheres, the core ex- 
pansion phase leads a significant part of the SIDM halo 
to a state close to that of a non-singular sphere. This 
is for example the case of the halo at maximum expan- 
sion in Fig. 1 (squares) where r c sa (9a 2 /AirGpc) 1 ^ 2 and 
r c is roughly equal to the radius where C is still posi- 
tive. Interestingly, the state of this halo is close to the 
non-equilibrium state of minimum entropy described by 
Lynden-Bell & Wood (1968). At the radius where the 
velocity dispersion (temperature) is still roughly constant 



(the radius of the isothermal sphere), the density has de- 
creased from the center by approximately a factor of 30 
and E and C are already negative. Systems like this are 
susceptible to the gravothermal catastrophe and this is 
what we observe it happens in our simulations. 

For our monolithic NFW halo, the maximum constant 
density core radius is reached approximately in 2.6 x 10 9 
yr. This time scale agrees with the estimated core relax- 
ation time, t TC , by KW for a <~ 0.8 in their dimensionlcss 
units or, equivalently, o~dm ~ 5.0 cm gr . As mentioned 
above, had we used V12 instead of the particle velocity, 
a higher scattering rate would have been expected. The 
factor 9/5, however, is lower than the factor 3 that KW 
claim they measure. The difference is probably due to the 
way in which we compute local density. If we estimate 
the core radius r c by the point where the density drops 
to 1/4 the central density, then r c <~ 0.4rNFW- Using the 
same definition for r c , finds r c ~ 0.6r p , where r p is the 
radius where the velocity dispersion peaks for a Hernquist 
density profile. Because the radius r p is close to tnfw for 
a NFW density profile, we would have expected, for pure 
analogy with Burkert's results, that r c ~ 0.6rNFW- We 
see from Figure 1, panel (b), that the peak is broad and 
is slightly lower than tnfw- It is not thus difficult to see 
why the core size is not as large as in the Hcrnquist's case. 
We let the simulation run for another 8.5 x 10 9 yr and 
compared the rate at which the core collapse is produced 
with that predicted by KW (see their Fig. 2). In the KW 
halo, for which t IC ~ 2.6 x 10 9 yr and a ~ 0.8, the core 
collapse develops at a faster rate than the one we observe 
in our halo. Our core density evolution is more akin to a 
KW a simulation with a value between 1.0 and 3.0. It is 
not clear what this difference is due to, but this difference 
could also be explained by the fact that we are using a 
NFW density profile instead of a Hernquist one. In any 
case, we find a good agreement between our results and 
Burkert's and KW results. 

2.2. Cosmological simulations of SIDM halos 

The set of simulations performed to analyze the struc- 
ture of halos in a SIDM cosmology uses the multiple-mass 
variant of ART (Klypin et al. 2001) in order to increase 
the mass and spatial resolution in few selected halos. All 
models assume a ACDM power spectrum with total matter 
density and cosmological constant, in units of the critical 
density, of fi m = 0.3 and Q\ = 0.7, a Hubble constant of 
h = 0.7, in units of 100 km s" 1 Mpc -1 , and a cr 8 = 1.0. 
Here erg represents the rms of mass fluctuations estimated 
with the t op-hat window of radius 8h _1 Mpc. 

We first perform two low-mass resolution simulations, 
one in a 100 h _1 Mpc box and the other in a 12.5 h _1 Mpc 
box. Both simulations were run with 64 3 particles in a 
grid initially consisting of 256 3 cubic cells. Two cluster- 
sized halos from the 100 h _1 Mpc box simulation and five 
MW-sizcd halos from the 12.5h _1 Mpc box simulation were 
selected for being rcsimulatcd with higher resolution. To 
increase resolution in the chosen halos, particles within 
2.5 virial radii from the halo centers were traced back 
to the initial epoch and new initial conditions were then 
generated by using the multiple-mass scheme described in 
Klypin et al. (2001). Three different mass levels with 1, 
8, and 64 times m p were used, where m p is the mass per 
particle at the finest mass level (see Table 1). Only parti- 
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Table 1 



Lbox 


m p 


Resolution 


^0 


a 




"max 


A' 


Ncol 


Halo name tag 


(h^Mpc) 


(h-^o) 


(h^kpc) 


(cm 2 gr~ r 


1 




(kms -1 ) 


(io- 2 ) 






100.0 


5.0 x 10 9 


3.0 


0.0 


0.0 


4.5 x 10 14 


1280 


4.43 


0.00 


CZlo.o 


100.0 


5.0 x 10 9 


3.0 


0.0 


0.0 


5.0 x 10 14 


1230 


1.80 


0.00 


C72o.o 


100.0 


5.0 x 10 9 


6.1 


0.1 


0.0 


4.4 x 10 14 


1285 


4.57 


0.60 


C71 .i 


100.0 


5.0 x 10 9 


6.1 


0.5 


0.0 


4.5 x 10 14 


1299 


4.88 


2.37 


CIIq.5 


100.0 


5.0 x 10 9 


6.1 


3.0 


0.0 


4.4 x 10 14 


1336 


4.62 


12.54 


Clh.o 


100.0 


5.0 x 10 9 


3.0 


3.0 


0.0 


3.6 x 10 14 


1548 


2.17 


30.03 




100.0 


5.0 x 10 9 


6.1 


0.5 


1.0 


4.4 x 10 14 


1276 


4.47 


0.35 


C^lo.5/V 


100.0 


5.0 x 10 9 


6.1 


1.0 


1.0 


4.4 x 10 14 


1270 


4.45 


0.64 


C'li.o/v 


12.5 


9.7 x 10 6 


0.2 


0.0 


0.0 


1.4 x 10 12 


215 


3.82 


0.00 


GIq.o 


12.5 


9.7 x 10 6 


0.2 


0.0 


0.0 


2.7 x 10 12 


237 


7.21 


0.00 


G2o.o 


12.5 


9.7 x 10 6 


0.2 


0.0 


0.0 


8.7 x 10 11 


167 


6.90 


0.00 


G3o.o 


12.5 


9.7 x 10 6 


0.4 


0.1 


0.0 


1.3 x 10 12 


214 


4.18 


0.33 


Gl .i 


12.5 


9.7 x 10 6 


0.4 


0.5 


0.0 


1.3 x 10 12 


217 


4.45 


1.32 


Glo.5 


12.5 


9.7 x 10 6 


0.4 


0.5 


0.0 


2.6 x 10 12 


240 


7.45 




Gr2o 5 


12.5 


9.7 x 10 6 


0.4 


0.5 


0.0 


1.1 x 10 12 


193 


1.85 




G4 . 5 


12.5 


9.7 x 10 6 


0.4 


0.5 


0.0 


1.1 x 10 12 


178 


5.53 




G5 .5 


12.5 


9.7 x 10 6 


0.4 


3.0 


0.0 


1.2 x 10 12 


245 


4.50 


9.41 


GI3.0 


12.5 


9.7 x 10 6 


0.1 


3.0 


0.0 


2.3 x 10 12 


273 


5.70 


15.03 


G2 3 . 


12.5 


9.7 x 10 6 


0.4 


3.0 


0.0 


8.3 x 10 11 


161 


6.02 


4.02 


G33.0 


12.5 


9.7 x 10 6 


0.4 


0.5 


1.0 


1.3 x 10 12 


212 


4.28 


0.77 


Gl .5/y 


12.5 


9.7 x 10 6 


0.4 


0.5 


1.0 


2.7 x 10 12 


239 


7.53 


0.95 


G2().5/V 


12.5 


9.7 x 10 6 


0.4 


1.0 


1.0 


1.3 x 10 12 


212 


4.31 


1.34 


Gli.o/y 



cles at this level were allowed to collide. Final halos have 
a few 10 5 particles each and contamination within their 
virial radii by heavy particles was kept below 1%. We 
notice that initial conditions do not change when obm is 
varied. 

The bound density maxima (BDM) group finding algo- 
rithm (Klypin et al. 1999) was used to identify halos and 
subhalos in all simulations. The BDM algorithm finds po- 
sitions of local maxima in the density field smoothed at the 
scale of interest and applies physically motivated criteria 
to test whether a group of particles is a gravitationally 
bound halo. In the BDM algorithm, there is a param- 
eter called the rejection velocity limit vi which monitors 
whether a particle is bound or not. For standard CDM 
simulations, in which a virializcd halo has a NFW like 
density profile, we usually set vi = v e , where the escape 
velocity v e is analytically evaluated in BDM assuming the 
halo follows a NFW density profile. We decided to turn 
off the rejection velocity switch because density profiles 
of SIDM halos are not supposed to be fitted by a NFW 
profile. For subhalos the switch is kept on for two reasons: 
first, a significant fraction of unbound particles is expected 
to populate these subhalos, and second, because these ha- 
los are composed of at most hundreds of particles, their 
profiles in most cases are well fitted by a NFW profile. 

2.3. Density Profiles 

A summary of the parameters of the simulations of ha- 
los are presented in Table 1. The mass of the particle in 
the high-resolution region m p is shown in column (2) and 



is only a funtion of the size of the box (col. [1]) because 
the cosmological model and the number of mass levels are 
fixed. Column (3) shows the formal force resolution, mea- 
sured by the size of a cell in the finest mesh. The col- 
lisionlcss CDM simulations have a higher resolution than 
the collisional ones because simulated halos in the former 
case reach a higher central density than in the latter one. 
Values of the parameters of the SI cross section, <7 and 
a, are presented in columns (4) and (5). The mass M vir 
within the virial radius, defined as the radius at which the 
average halo density is S times the background density ac- 
cording to the spherical top-hat model, is shown in column 

(6) . 5 is a number that depends on epoch and cosmologi- 
cal parameters (O to ,Oa); for a flat ACDM model, 6 ~ 337 
at z = 0, where z is the redshift. The maximum circular 
velocity defined as 

( GM{<r) \ 1/2 

»max = I I , (3) 

V / max 

where G is the gravitational constant and M(< r) is the 
mass contained within the radius r, is placed in column 

(7) . The spin parameter A' as defined by Bullock et al. 
(2001) is shown in column (8) (see §3.4). In column (9) 
N co i denotes the total number of collisions until z = 
divided by the total number of particles within H V i r ; as 
mentioned above, only particles at the finest mass-level of 
refinement arc allowed to collide. This number could not 
be obtained for those halos (horizontal dots) that were 
simulated simultaneously in only one run because the col- 




Fig. 3. — Density profiles of the fiducial cluster- (solid lines) and MW-sized (dashed lines) halos simulated with different values 
of <jdm (see top of each panel). V is the V12 in unities of 100 kms -1 . Dotted curves are the profiles of the CDM case (panel a) 
plotted for comparison. Dashed curves in panel (a) are the best NFW fit to the profiles shown there. Radius is given in unities of 
the virial radius R vir . 



lision counting refers to the whole region selected to be 
rcsimulated at high resolution. In column (10) we give the 
halo name tag for each one of the simulated halos. La- 
bels Gin denote cluster-sized halos (n goes from 1 to 2), 
whereas labels Gn denote MW-sized halos (n goes from 1 
to 5). The numbers n = 1, 2 for halos CI, and n = 1, 2, 3 
for halos G indicate different simulations aimed to have 
halos with different mass assembly histories. MW-sized 
halos numbers 4 and 5 turned to appear in the same sim- 
ulation as the halo G2, and we have included them for 
completness. Both, the cluster- and MW-sized halos have 
been simulated for a variety of (ao,a) values; the lower in- 
dex in the name tag indicates the used combination. Note 



that only halos with n = 1 were simulated with all the six 
{<Jo,a) pairs of values; they are our fiducial halos. 

3. HALO STRUCTURE 

The six panels of Figure 3 show the density profiles of 
the fiducial cluster-sized halo (Cll, solid lines) and MW- 
sized halo (Gl, long-dashed lines) for the different {(Jo, a) 
pair values considered here. The fiducial halos Gl and 
Cll were selected in the low-resolution CDM simulations 
with concentrations typical for their masses in an attempt 
to simulate halos with typical MAHs, since it is known 
that the concentration of the halo is related to its MAH 
(Avila-Reese et al. 1998, 1999; Wechsler et al. 2001). Un- 
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fortunately, in the high-resolution CDM simulation, the 
halos turned out to be more concentrated (cnfw = 10 m ~ 
stead of 6, and cnfw = 17.2 instead of 13, respectively). 
Besides, a major merger in the halo Cll that occurs at 
a < 1 in the low-resolution CDM simulation is delayed to 
a ;> 1 in the corresponding high-resolution simulation; this 
will shift the core collapse phase to an epoch a > 1 when 
(7m = 3.0 cm 2 gr _1 , as we will discuss below. In Figure 3, 
for reference, the density profiles of the two CDM models 
shown in panel (a) are also plotted with dotted lines in the 
other panels. The best NFW fits to the CDM profiles are 
shown with short-dashed lines in panel (a). 

From a first inspection of Figure 3, we conclude that 
significant differences in the density profiles (flattening) 
among the halos with different values of (<to, a) apply only 
to those regions with radii smaller than 0.03 — 0.05 R„i r . 
The exception are the halos with (a ,a) = (3.0,0.0); in 
this case SI affects dramatically the inner structure at radii 
even larger than 0.05R„i r . The parameter which defines 
the evolution under SI is the number of collisions per par- 
ticle after a given time t, 

N col oc p {a dm | V12 1) t oc p<7 DM cr 3D t, (4) 

where u^d is the 3D velocity dispersion. The maximum 
03D is proportional to v max . From all of our runs we see 
that if N co i at the present epoch is smaller than « 2 — 5 
(see Table 1), then the halo is either in the core expansion 
phase or has just had the gravothermal catastrophe trig- 
gered, being the core shrinkage still negligible. When o~dm 
is constant, the cores of low velocity halos are on average 
less influenced by SI than the high velocity ones. On the 
other hand, when a dm °c I/V12, we expect halos of differ- 
ent sizes to have roughly similar core densities. These pre- 
dictions are closely obeyed as seen in panels (b), (c), and 
panels (e) and (f) of Figure 3, respectively, or in Figure 4 
below. Nevertheless, this simple reasoning applies strictly 
only for monolithic halos. The cosmological merger history 
may dramatically influences the z = structure of halos 
with significant SI. This is the case of the cluster-sized 
halo in the simulation with gdm = 3.0 cm 2 gr _1 , CZI3.0 
(panel d). This halo should have evolved faster under the 
influence of SI than the MW-sized halo Gl3.o- However, 
this did not happen; while halo GI3.0 is in an advanced 
stage of core collapse, halo C71 3 . is just entering into this 
phase. The core collapse in the latter case seems to have 
been delayed by a recent major merger. We have simu- 
lated a second cluster-sized halo with <jdm= 3.0 cm 2 gr _1 , 
C72 3 . , which does not suffer late major mergers. The core 
collapse phase at a = 1 in this halo is well advanced and 
its central density and concentration are higher than those 
of the corresponding MW-sized halo (see Fig. 4 below). 

An important effect to note is that the core evolves fast 
after it reaches its maximum expansion. The central den- 
sity as well as N co i (see eq. 4) increase quickly. This is 
why N co i at z = for the halo GI3.0 is similar to that 
of the halo C7l3.o(see Table 1) despite that the latter has 
a v max 5.5 times larger than the former. We recall that 
halo C7l3.odelayed its core collapse due to a late major 
merger. Owing to the transient nature of the SIDM halo 
cores, a comparison of core parameters at a given epoch 
(e.g., z — 0) for different values of odm and different halos 

2 Aside from ctq we also need to know the concentration and the total 



masses and MAHs may lead to misleading conclusions. 

In the following, we analyze in more detail the density 
profiles of the cluster and MW-sized SIDM halos in the 
light of previous works and then explore the influence of 
the MAH on their evolution. 

3.0.1. Cluster-sized halos 

As mentioned above, it turned out that the chosen 
cluster-sized halo, which in the low-resolution CDM sim- 
ulation was relaxed at a = 1.0 (a = 1/(1 + z)), is in the 
high-resolution simulation about to suffer a major merger. 
The time when this merger happens seems to be a function 
of a dm- We find that the halo Cll with {a ,a) = (3.0,0.0) 
suffers this merger before it reaches the present time, while 
Cll halos with other <jdm values suffer this merger soon 
after a = 1.0. This is why in Figure 3 we compare the den- 
sity profiles of halos Cll at time a = 1.1, an epoch when 
the Cll halos with different values of a dm have already 
experienced their last major merger. 

We find the same general trend observed by Yoshida et 
al. (2000b) at z = 0: as odm increases, the core radius, 
r c ,_i, defined in this paper as the radius where density 
profile slope becomes steeper than —1, increases and the 
central density, p c ,—i, decreases. However, differences be- 
tween Yoshida ct al. results and ours do exist: in our 
simulations there is a (Jdm (for the constant odm case) 
for which the trend is reversed. For example, r Ci _i for halo 
CZI3.0 is not greater than the corresponding r Ci _i of halo 
C71o.5, and for halo C723.o, r c ,-i is even much smaller. 
The trend of the core radius or central density with udm 
at z — is reversed owing to the core collapse. The halo 
with (Jdm = 10.0 cm 2 gr _1 of Yoshida et al. is also un- 
dergoing a core collapse; however, at z = the central 
density of this halo is still larger than the central den- 
sity of their halo simulated with <tdm — 1-0 cm 2 gr _1 . Is 
the rate at which the core shrinks and the central density 
grows up in the stage of core collapse slower in their case 
than in ours? We have measured this rate in both their 
cluster and our cluster halo GZI3.0 and have found that 
the rates are actually similar. The difference between ours 
and Yoshida's et al. results may be found in the minimum 
central density reached by the halos. In their case, this 
density is lower than in our case by a factor of three. This 
is probably why the soft core at z = of their cluster halo 
simulated with odm = 10.0 cm 2 gr _1 have not become yet 
smaller than the corresponding one to the simulation with 
o dm = 1-0 cm 2 gr _1 , in contrast to our results where for 
o dm = 3.0 cm 2 gr _1 the soft core is already smaller than 
for o dm = 0.5 cm 2 gr _1 . 

Since halo CII3.0 assembles most of its mass through 
late major mergers, the triggering of the core collapse 
phase, which is almost inevitably present in almost all ha- 
los with (cr , a) = (3.0,0.0), is delayed, and at a = 1.1 still 
has a prominent soft core, although it is already shrinking 
because of the core collapse. After the last major merger 
took place, little mass is accreted by the halo in the next 
3.5 x 10 9 yr; it is thus tempting to compare the core evolu- 
tion since that epoch with the predictions of the monolithic 
case of KW. The halo C71 3 . would have a a that lies in the 
interval [3.5,4.0] 2 and a dynamical time td yn , as defined 
by KW, of 1.7 x 10 9 yr. From a = 0.85, which is the epoch 

mass in order to compute a. We have taken these parameters from 
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where the maximum core expansion occurs, to a = 1.10 
we see that the core shrinks by a factor of 1.3. KW predict 
for this model a shrinking factor of about 1.4-1.5. It is 
not clear how far we can take the comparison of the SIDM 
cosmological scenario with the monolithic case, but here 
we have a case in which consistent results are obtained. 
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Fig. 4. — Central density, core radius and core fraction vs. 
v max for cluster- and MW-sized halos simulated here for dif- 
ferent values of (7dm (see the corresponding symbols in the 
middle panel; udm is in units of cm 2 gr _1 and V is V12 in units 
of 100 kms^ 1 ). Central density p c ,-i and core radius r c> _i are 
deffined as the density and radius where the profile slope be- 
comes steeper than — 1, respectively. Dashed and dotted lines 
in the two upper panels are the average and the 2 — a disper- 
sion inferred from observations of dwarf and LSB galaxies and 
cluster of galaxies (see §5.2 for details). 



3.0.2. Milky Way-sized halo Gl 

The fiducial MW-sized halo Gl resulted more concen- 
trated (see above) and with a MAH more biased to an early 
halo-assembly than the average corresponding to its mass. 
Moreover, the late MAH of this halo is nearly smooth, 
without violent major mergers since a — 0.5. 

The non-monothonic trend of the core radius or the cen- 
tral density with a dm at z — sliglthy outlined by the 



fiducial cluster-sized halo is clearly present in the fiducial 
MW-sized halo (Fig. 3). This is much better appreci- 
ated in Figure 4. Here we plot v max versus central density, 
p c ,-i, measured at the radius f c ,-i where the profile slope 
becomes steeper than —1 (a), the core radius, r Ci _i (b), 
and the core radius fraction, r C) _i/R„j r , (c). Note that 
halos with (<7n,a) = (0.1,0.0) and (0.5,0.0) tend to lie in 
the same region; i.e., an increase in a docs not necessar- 
ily produce an increase in core radius at z = 0. In the 
framework of a monolothic halo this could be explained 
as follows: because the core relaxation time for halos with 
((To, a ) = (0.1,0.0) is close or greater than the Hubble time, 
we still find these halos in their core expansion phase at the 
present time, whereas halos with (cto, a) = (0.5,0.0) are at 
the onset of the core collapse phase. These facts conspire 
to give similar core parameters in both cases. Halo GI3.0 
on the other hand, is well inside the core collapse phase 
and this is why r Ci _i is small. It is interesting to note 
that SI introduces a peculiar epoch in the history of the 
universe at which the dark halos change their inner struc- 
ture from an extended smooth core to a very cuspy center. 
For the case a dm oc 1/vi 2 , this epoch does not depend 
too much on the halo mass or v max , and it is determined 
mainly by the value of <jq. However, violent mergers and 
other effects may delay significantly the core collapse phase 
for some halos and subhalos (see §3.1.3). 

In Figures 5 and 6 we plot the evolution of the density 
and 3D velocity dispersion profiles, and of the heat capac- 
ity profile, respectively, of the core collapsing halo GI3.0 
from a = 0.4 to a = 1.0 (left panels). Proper units were 
used. As in the case of halo C71 3 . , one could ask whether 
the strong evolution of p c is expected or not. To address 
this question we first notice that only 20% of the present 
total mass is accreted since the last major merger, which 
occurs at a ~ 0.5. Therefore, treating the evolution of this 
halo from a — 0.5 to a — 1 as an isolated case seems to be 
a fair approximation 3 . As we did for the fiducial cluster- 
sized halo, we compute the corresponding a and dynamical 
time; they are 1.9 and 6.4 x 10 8 yr, respectively. From the 
KW simulations of monolithic halos, this model changes 
its core radius by a factor of four in approximately 5-6tdy„; 
this same change is seen in halo GI3.0 but in about lOtd^. 
We can explain this difference in the core evolution rate 
as being due to the infalling and merging material in the 
cosmological halo. This also could explain why the in- 
ner region where the heat capacity C is positive increases 
during the expansion phase in halo GI3.0 (Fig- 6) unlike 
what happens with the monolithic halo in §2.1 (Fig. 1). 
The subsequent core collapse phase is characterized by a 
reduction of the region where G > 0. 

The core collapse of halo GI3.0 (and G23.0, see below) 
is at odds with the results of Dave et al. (2001). They 
ran a simulation of 128 3 particles in a 4 h _1 Mpc box for 
the same cosmological parameters used here. Their power 
spectrum is normalized to eg = 0.8. The most massive 
halo formed in their simulation has a mass of 6.2 x 10 11 Mq 
and consists of about 1.7 x 10 5 particles. None of the ha- 
los in the two SI models that they consider, (0.56,0.0) and 
(5.6,0.0), shows signs of being in the core collapse phase. 



its correspondent CDM cluster halo. This is a good approximation because once the cluster forms, which we can establish at the last major 
merger, the parameter tnfw does not change much. 

3 We later realized that what matters is not only the amount of mass but how this mass is accreted. Smooth infall of material to halo 
outskirts does not seem to alter significantly the dynamical core SI evolution 
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Fig. 5. — Evolution of the density and the 3D velocity dispersion profiles of the Gl (left) and G3 (right) halos in proper units 
and for (<tdm, a) — (3.0, 0.0). In both halos the core reaches its maximum size a little later after the major marger occurs, around 
a = 0.5 for GI3.0 and a ~ 1.0 for G33.0; for times earlier than this, the central densities of the major progenitors of the halos 
are higher. The shape of the density profiles of their corresponding CDM halos do not change after the major mergers occur, the 
subsequent infall of material serves essentially to increase the proper size of the halo. Because a very small amount of material 
reach the core (see also the paper by Zhao et al. 2002), we can safely say that the inner regions of the SI halos evolved as if they 
were isolated; in other words, the number of collisions inside the core overwhelms the number of particles that fall into it. So, it 
is not very surprising, after all, that we also see a core collapse in MW-sized cosmological halos. 



On the contrary, the halos have extended shallow cores 
which increase as <jdm increases. Dave et al. suggested 
that the accretion of dynamically hot material, intrinsic 
to the CDM cosmological simulations, is so important to 
the halo core evolution that predictions for SIDM from 
monolithic halos does not apply at all to the cosmological 
setting. 

3.0.3. Halos with different assembling histories 

In the following, we test the posibility that our differ- 
ences with Dave et al. (2001) might reside in the halo 
MAHs. It is well known that CDM halos form from a va- 
riety of MAHs which give rise to a variety of present-day 
structures (e.g., Avila-Reese et al. 1998,1999). As men- 
tioned above, the fiducial halo assembled early and ended 
highly concentrated. We need to simulate now a halo with 
a low concentration and a late halo-assembly. We chose 
two MW-sized halos with these characteristics, G2o.o and 
G3o.o, from the same low-mass resolution simulation as the 
one from which the halo Glo.o was taken; G23.0 and G33.0 
are the corresponding SI halos with ((T , a) = (3.0,0.0). As 
for the fiducial halo, in the high resolution simulation the 
halos resulted with a higher concentration and a less ex- 
tended MAH than in the low resolution one. In any case, 
both halos are less concentrated and have a more extended 
MAH than the fiducial one. 

To our surprise, halo G2 3 . also undergoes a fast core 



collapse, no matter that this halo accretes a little more 
than half of its present mass from a = 0.5 to a = 1.0! An 
analysis of the MAH of the corresponding CDM halo shows 
that the last major merger occurs as early as a ~ 0.3; after 
this epoch the mass is accreted smoothly and/or in small 
chunks. A recent paper by Zhao et al. (2002) shows that 
the build-up of CDM halos can be divided in time in two 
phases: an early one characterized by violent mergers in 
which the inner halo density profile is significantly altered, 
and a late one characterized by a smooth acretion of mate- 
rial in which the accretted material retains its low binding 
energy and is added to the outer part of the halo. De- 
pending on the MAH, a halo can go from one phase to the 
other early or late in time. Halos Glo.o and G2o.o appear 
to have typical transition epochs, although the latter ac- 
cretes (smoothly) much more mass at late epochs than the 
former one. On the other hand, halo G3o.o seems to have 
the transition at z ~ 0. This halo assembles most of its 
mass by recent major mergers. 

Besides of the direct dynamical effect of the mergers on 
the structure of halo G33.0, its core collapse may be de- 
layed by conduction effects. The major mergers heat the 
halo periphery, producing an increment of the outer veloc- 
ity dispersion (temperature), as it seen in the right panel 
of Fig. 5 (dotted line). In this case, the halo at a = 0.4 was 
already starting its core collapse phase. The heat capacity 
also increases due to the major mergers (Fig. 6). Then, 
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Fig. 6. — Evolution of the heat capacity, C = dU/dT, profiles for halos GI3.0 (left panel) and G33.0 (right panel). Notice how 
the size of the region where heat capacity is positive increases in both halos as the cores of the halos approach their maximum 
radius. As the core collapse proceeds in halo GI3.0 the region where C is positive shrinks. This is different to what we see in the 
subsequent evolution of C in halo GI3.0, where the core collapse is absent; in this case the size of the G > region even increases 
as time passes. 



heat starts flowing into the C > core and produce its ex- 
pansion (see §2.1). This process may last a significant pe- 
riod of time, until the heat conduction isothcrmalise most 
of the halo. We let the simulation of halo G33.0 run for 
another Hubble time and to our surprise even at a — 1.6, 
about 20 Gyr after the Big Bang, the halo has not under- 
gone a core collapse and the heat capacity increased even 
slighlty. We should note that this halo is also peculiar be- 
cause close to it there is a galaxy-sized more massive halo. 
We notice, on the other hand, that because halo G3 3 . is 
smaller (v max = 163 kms^ 1 ) than halos GI3.0 and G23.0, 
we expect the SI evolution of this halo to be slower than 
that of halos GI3.0 and G23.o- 

Our results show that, for both the cluster- and MW- 
sized halos with {?dm , o) — (3.0, 0.0), the core collapse is 
delayed until the stage of major mergers ends. We further 
simulated the halo Gl with a still higher o~r>M (= 10.0 
cm 2 gr -1 ) and check, in fact, that the core collapse begins 
only after the last major merger. Therefore, we conclude 
that the dynamics of SIDM halos with a high cross sec- 
tion is strongly affected by violent merging processes rather 
than by smooth mass aggregation. A continuous injection 
of "dynamical heat" , able to either heat and expand di- 
rectly the core or rise the temperature of the halo periph- 
ery to above the core's temperature, seems to be possible 
only when the halo grows by violent mergers. In summary, 
when SI is efficient {o~dm = 3.0 cm gr , for example) the 
inner structure of halos may depend strongly on the halo 



merging history. The onset of the core collapse phase oc- 
curs immediately after the core relaxes, in agreement with 
the finding of KW (see also Quinlann 1996; Burkert 2000), 
but this process can be delayed by subsequent major merg- 
ers as we actually see in some of our halos. Can this explain 
our differences with Dave et al. (2001) results? 

3.1. Concentrations 

A more quantitative way to measure differences in the 
structure among halos is using concentration parame- 
ters. In Figure 7 we plot the C1/5 and cm concentration 
parameters 4 versus virial mass M vir for halos (large sym- 
bols) and subhalos with more than 1000 particles (small 
symbols) from our different simulations. To avoid too 
much overlapping, we divided Figure 7 in two sequences of 
panels. Halos with (a ,a) = (0.0,0.0), (0.1,0.0), (3.0,0.0), 
and (1.0,1.0) are shown in left panels while halos with 
(o-Q,a) = (0.5, 0.0) and (0.5,1.0) are shown in right panels. 
We also plot in Figure 7 the density measured at 0.03R„i r 
(Pc.03) versus M„j r . Because our high-resolution simula- 
tions were designed to be focused on a specific (cluster- or 
MW-sized) halo, most subhalos lie within the virial radii 
of these halos (satellites) . The mass of subhalos is defined 
as the minimum between the virial mass and the truncated 
mass (see Avila- Reese et al. 1999). 

From Figure 7 one sees that halos with different (cr , a) 
values have similar C]y 5 parameters, except for halos GI3.0 
and G2 3 . , in which the region where l/5th of M vir is con- 

4 c 1 / 5 is defined as the ratio between the virial radius K v i r and the radius where 1/5 of the total halo mass is contained (Avila- Reese et al. 
1999). cm is defined as 27 times the ratio between the mass at ri n = 8.5 kpc 

"max 

/220kms- 1 and the mass at r ou t = 3xr; n (Dave et al. 2001). 
These definitions of the concentration parameter are independent from the particular fitting applied to the halo density profile. 
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tained is greatly affected by the core collapse, and halo 
G33.0, whose peculiar MAH makes it keep a relatively 
long-lived shallow core (see §3.1.3). A similar behavior is 
observed for p Cj .o3- The parameter Cm reflects with better 
accuracy the inner matter distribution of halos. Note that 
the Cm parameter of subhalos with more than 1000 parti- 
cles from the MW-sized (<7o, a) — (0.5, 0.0) simulations are 
similar to those from the (cro,a) = (0.56,0.0) simulation 
of Dave et al. (2001). Halos with higher have smaller 
cm concentrations, for a given a. Again, this trend can be 
reversed at the high oo end because of the core collapse. 
Subhalos follow roughly the same behavior of halos in all 
the three parameters of Figure 7 as a function of mass. 
However, before drawing any conclusion we should bear 
in mind that (i) subhalos are resolved only with hundreds 
of particles and (ii) that they can result strongly affected 
by the self-interacting environment of the halo where they 
are embedded. 

The lack of mass resolution does not allow us to ex- 
plore in detail the effect of SI on the inner dynamics and 
structure of subhalos. In order to evaluate this problem, 
we analyze the MW-sized halos obtained in a low reso- 
lution (64 3 particles) simulation with (ao,a) — (3.0,0.0). 
This is the same simulation from where MW-sized halos 
were taken. The SIDM MW-sized halos in this simula- 
tion have only few thousands of particles and we find that 
their inner structures are not at all similar to those of 
their high-resolution counterparts. Therefore, we can not 
be conclusive about the inner structures of the subhalos 
in the high-resolution simulations presented here because 
they have at the most only a few thousands of particles. 
Nevertheless, our simulations allow us to explore the ef- 
fect the hot particles of the host halo has on the overall 
structure of subhalos. The C1/5, Cm, and p c ,.03 parame- 
ters of the subhalos from the cluster-sized halo simulation 
CZI3.0 are smaller than those same parameters of the iso- 
lated MW-sized halos from the low resolution simulation. 
In both cases the halo masses and number of particles are 
similar. This difference seems to have only one explana- 
tion: the structure of subhalos (satellites) is affected by 
the environment of the hot host halo. Subhalo particles 
are being ejected continuously because high speed particles 
from the hot host halo collide with the cooler subhalo par- 
ticles, exchanging orbital energy by internal energy; in the 
extreme case, the subhalo can be evaporated completely 
(Gnedin & Ostriker 2002). The survivor subhalos are in- 
deed "puffy" with density profiles much flatter than the 
isolated halos. 

More evidence about the influence of SIDM halos on 
the structure of their subhalos comes from comparing 
subhalos of simulation Gli.o/y with subhalos of simula- 
tion GI3.0 (Fig. 7). Subhalos from the former simulation 
(squares) have an effective a dm of k 3 — 5 cm 2 gr _1 , 
i.e., similar to the udm of the latter simulation (lowest 
starred symbols); however, they are located in diferent 
regions in Figure 7 (the same happens for subhalos with 
less than 1000 particles). The difference between both 
simulations is in the effective a dm value of the host 
halos: ~ 0.5 cm 2 gr _1 versus 3 cm gr , respectively. 
The evaporation effect of the hot SIDM halo GI3.0 is 
stronger than the one of the halo Gl 10 /y; in fact, in 
the former case most subhalos have been destroyed by 
the present time (see §4). A similar effect is observed 



in subhalos from the (<Jo,a) — (0.5,1.0) and (0.5,0.0) 
simulations. For the MW-sized halos of the former sim- 
ulation, the effective value of udm is ~ 0.2 cm 2 gr _1 
while for their subhalos it is odm ~ 1-5 — 3.0 cm 2 gr _1 ; 
i.e., 3-6 times larger than for subhalos of the latter sim- 
ulation. Despite the relatively high value of a dm of 
subhalos in the (0.5,1.0) model, their structures look 
similar to their CDM counterparts (low-mass resolution 
does not allow SI to work properly). However, subha- 
los from the simulation with a dm — 0.5 cm 2 gr _1 are 
"puffier" than those subhalos from CDM (right panels, 
Fig. 7). This can also be explained by an action of 
the hot halo environment on the dynamics of subhalos. 
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Fig. 7. — Concentration parameters Ci/ 5 and cm, and den- 
sity p c ,.03 at 0.03R„i r versus halo mass for halos (large symbols) 
and guest halos with more than 1000 particles (small symbols) 
from simulations with different values of a dm (see the cor- 
responding symbols in upper panels; same unities as in Fig. 
3). Results were divided in two sequences of panels in order 
to avoid too much overlapping. Tiny starred symbols in the 
left panels are for the low-mass resolution simulation in the 
12.5 h^Mpc box with (a ,a) = (3.0,0.0). Their structural 
properties are different to that of their high-mass resolution 
counterparts (bright starred symbols). Dashed line in the up- 
per panels is the linear regression found in Avila-Reese et al. 
(1999) for isolated ACDM halos in a 60 h^Mpc box. Dashed 
line in medium and lower panels are estimates of the cm — M. v i r 
and p c ,.03 — M„; r relations for isolated ACDM halos from the 
the same simulation of Avila-Reese et al. (1999). ca^ was 
calculated from cnpw using Fig. 6 in Dave et al. (2001). 



3.2. Ellipticities 

Another prediction by the SIDM cosmology is that the 
halo inner region, where a shallow core forms, should be 
close to spherical because collisions tend to reduce the de- 
gree of anisotropy. 

We have measured the ellipticities of halos using the 
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Fig. 8. — Ellipticities el and e2 as a function of radius are shown in left (halo Gl) and right (halo Cll) panels. As expected, SI 
halo are rounder than the corresponding CDM halos, being the effect greater for the cluster-sized halo. 



tensor of inertia. This is defined as 



Iij — XiXj /r 



(5) 



where the sum is over all particles within r V i r , xi (i = 
1, 2, 3) are the particle coordinates with respect to the halo 
center of mass, and r is the distance of the particle to the 
halo center. The ellipticities are then given by 



ei 



1 



Ai 

As' 



e2 



1 



A2 

As' 



(6) 



where A3 > A2 > Ai are the eigenvalues of the tensor 
of inertia. We evaluate the triaxiality parameter using 
the following formula (e.g.,Franx, Illingworth, & de Zecuw 
1990) 



\ 2 

A 3 



X 2 

A 2 



A3 _ A l 



(7) 



A halo is prolate (oblate) if T = 1.0 (T = 0.0). 

The ellipticities e\ and e2 as a function of radius are 
shown in Figure 8 for our fiducial MW- and cluster-sized 
halos for three models a dm = 0.0, 0.5, and 1.0/V r cm 2 gr _1 . 
We have denoted with solid and dashed lines el and e2, 
respectively. We see, as expected, that the differences in 
the ellipticities between SI and CDM halos are found in 
the inner parts. At large radii, the difference in el or e2 
values among the plotted MW-sized halos is lower than 
the difference between the values of the ellipticities el and 
e2. The opposite occurs for the cluster-sized halo. Note, 
moreover, that unlike Gl halos, SI Cll halos are rounder 
at all radii. 

Cores of SI halos are indeed rounder than their CDM 
counterparts, but they are not spherical. Unfortunately, 
ellipticities in real galaxies are ill determined, so it is not 



clear how much we can constrain the SI parameters with 
this kind of observations. Miralda-Escude (2002) claims, 
however, to have constrained the constant cross section to 
be less than ~ 0.02cm 2 gr _1 with the determination of the 
ellipticity of cluster MS 2137-23 through a lensing study. 
The question remains whether the ellipticity he derives for 
MS 2137-23 does disagree with the one found in a typical 
cluster-sized simulated halo with a dm > 0.02cm 2 gr~ 1 . We 
let this analysis for a subsequent paper. 

3.3. Angular Momentum Distribution 

For each halo we compute the total angular momentum 



as 



J = rrij rj x Vj 



(8) 



i=i 



where rj and Vi are the position and velocity of the zth 
particle with respect to the halo center of mass. We follow 
Bullock et al. (2001) and define a modified spin parameter 
A' to characterize the global angular momentum of a halo 



A' 



3d 



V2M vir V mr R vi . 



(9) 



where J V i r is the angular momentum inside the virial ra- 
dius R V ir, and V V i r is the circular velocity at radius R V i r . 
In Table 1, column (8), we show A' for all of our halos. No- 
tice that A' increases with a dm m both cluster- and MW- 
sized halos for constant a dm, from <jq = 0.0 to 0.5. Nev- 
ertheless, this effect is small, amounting to 10%, 16% and 
3% for halos Cll, Gl and G2, respectively. It measures 
roughly the difference in the number of particles with low 
angular momentum in inner regions between CDM and SI 
halos. Moreover, for those halos for which odm oc I/V12, 
the effective a dm is so small that a factor of two of change 
in cr passes unnoticed by A'. 
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max J' Jmax 

Fig. 9. — Specific angular momentum j distribution for halos Glo.o (circles) and GI0.5 (squares) panel (a) and for halos Cllo.o 
(circles) and Clio. 5 (squares) panel (b). Only cells with positive j were considered and the total mass M t was redefined to be the 
mass contained in all those cells with positive j. The line is taken from the analytical fit of Bullock et al. (2001) m(< j) = fij / 1 + j, 
where m(< j) = M(< j)/M vir and j = j/jmax, with /i = 1.25. Notice that differences between the CDM and SI m(< j) profiles 
are only at the low-j end. We see that our fiducial MW-sized halo is not well represented by the formula of Bullock et al., which 
agree with the findings by Chen & Jing (2002) . 



These results already suggest that significant differences 
in the angular momentum distributions (AMDs) of halos 
with and without SI are not expected. To compute the 
AMD we find first, for each particle, the specific angu- 
lar momentum component along J v i r , calling it j. We 
then divide the sphere of radius R V i r in spherical shells 
and each of these shells is in turn divided in four quad- 
rants. Differences in the number of particles between cells 
for most of them are below a factor of two. Only a small 
fraction of these cells have negative j values and we re- 
ject them in the further calculation of the AMD; notice, 
however, that these particles with negative j can have a 
significant effect on the angular momentum problem and 
the formation of bulges (van den Bosch et al. 2002). Cells 
are ranked according to their j value and M(< j) profiles 
are then built by counting the cumulative mass in cells 
with angular momentum smaller than j. Figure 9 shows 
the M(< j) profiles for halos Gl and Cll (panels a and 
b, respectively, where circles are for CDM and squares for 
o~dm = 0.5). Mass and angular momentum are given in 
M V i r and jmax units, respectively, where jmax is the max- 
imum value reached by j. We have also plotted in Figure 
9 the analytical form of M(< j) proposed by Bullock et 
al. (2001) with /j, = 1.25; the line intends to guide the 
eye, it is not a fit to our numerical results. The effect of 
the dynamical expansion of the core is seen only at low 
j. There is less mass, for a given j, in SI halos or, for 
the same fraction of mass, CDM halos have less rotational 
support. However, the difference between the CDM and 
the SI profiles is so small that SIDM can not account for 
the differences between the predicted angular momentum 



distribution of the CDM model and that determined for 
dwarf galaxies (van den Bosch, Burkert, & Swaters 2001). 

4. SUBHALO POPULATION 

In previous works on SIDM it was suggested that the 
number of subhalos in MW-sized halos should be lower 
than that predicted in CDM because subhalos are sup- 
posed to form with a lower concentration and/or because, 
the hot environment of a SIDM halo evaporates its subha- 
los. What do we see in our simulations? For halos GI3.0, 
G23.0, and G33.0 the evaporation effect is strong and most 
substructure is erased (see, for example, circles in Fig. 10, 
left panel). For o\dm=0.5 and 0.1cm 2 gr _1 (halos Glo.r 
and GI0.5), the effect is less pronounced but the number 
of subhalos is still less than in the CDM MW-sized ha- 
los, although at the low velocity end (v max = 20 kms -1 ), 
the number of subhalos becomes larger. A similar be- 
havior is seen in the halo Gl .5/y, where the effective 
o~dm of the MW-sized halo is s=s 0.2cm 2 gr _1 . The sub- 
halo v max — function of the halo Gli.o/y lies slightly above 
of the CDM one. The v max — function of subhalos seems to 
be sensitive to the effective a dm of the host halo rather 
than that of the subhalos. In any case, halos Glo.r, GI0.5, 
Gl .5/V-, and Gli.o/y have at least as many subhalos with 
Vmax ^ 25 kms -1 as the CDM MW-sized halo. Thus, we 
can conclude that the substructure problem (Klypin et al. 
1999; Moore et al. 1999) is not solved by a SIDM model 
with a < 1 and reasonable values of <7o, at least at the low 
Vmax end of the velocity function. 

The explanation of why the survival time of subhalos 
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in a SIDM scenario is comparable or even longer than the 
one in the collisionlcss CDM cosmology resides probably 
in the fact that the host halo tidal force, responsible for 
the disruption of substructure along with dynamical fric- 
tion, in the center is not as strong as in the CDM case. 
There may also intervene a numerical effect: the cores of 
subhalos, where the action of the self-interaction is impor- 
tant, are not resolved in the simulations in such a way that 
subhalo density profiles are very similar to those of CDM 
ones. Notice, however, that this argument applies only 
to relatively small effective <jdm values: if <jdm is very 
high, subhalo concentrations can be modified significantly 
by particle evaporation. 

For completeness we have also plotted the cumulative 
v max — function for the fiducial cluster-sized halo Cll with 
the different values of (a ,a) (Fig. 10, right panel). The 
total number of subhalos within the SIDM cluster-sized 
halos (a sphere of 1.5 h _1 Mpc radius was used) is only 
slightly smaller than for the CDM halo, excepting the cases 
with (aQ,a)= (3.0, 0.0) and (0.5, 0.0). It is interesting 
to notice that the evaporation constraint on the Funda- 
mental Plane of elliptical galaxies within the clusters set 
namely a maximum value of <jdm ~ 0-5 cm 2 gr _1 (when 
a = 0) according to Gnedin & Ostriker (2001). Regarding 
the shape of the v max — function, in the range of 200-300 
kms -1 , the SI cluster-sized halos have more subhalos than 
the corresponding CDM halo but the contrary happens at 
the low-v max end, SI cluster-sized halos have less subhalos 
than the CDM halo. This could be attributed to a longer 
survival time of massive subhalos. We do not expect a 
serious observational conflict related to the subhalo popu- 
lation in SIDM cluster-sized halos with (ao,a)= (0.1,0.0), 
(0.5,1.0) and (1.0,1.0). 

5. VIABILITY OF THE SIDM COSMOGONY 

Spergel & Steinhard (2000) proposed the SIDM model 
in an attempt to overcome the apparent discrepancy be- 
tween the cuspy halos predicted by CDM and the shallow 
halos inferred from observations of dwarf and LSB galax- 
ies (Moore 1994; Flores & Primack 1994; Burkert 1995; 
de Blok & McGaugh 1997). This discrepancy has been 
challenged by Swaters, Madore, & Trewhella (2000), van 
den Bosch, Robertson, Dalcanton, & de Blok (2000), and 
van den Bosch & Swaters (2001) who concluded that the 
spatial resolution of observed rotation curves was not suffi- 
cient to put constraints on cosmological models. However, 
more recent observational studies, using the highest sensi- 
tivity and spatial resolution in HI, Ha and CO lines show 
that halos of dwarf and LSB galaxies are actually less con- 
centrated in the center than the prediction of the CDM 
model (Cote, Carignan, & Freeman 2000; Blais-Oucllette, 
Amram, & Carignan 2001; Bolatto et al. 2002; Amram & 
Garrido 2002; de Blok ct al. 2001a, 2001b; Marchesini ct 
al. 2002). Also, there are pieces of evidence of soft cores in 
normal disk galaxies (e.g., Corsini et al. 1999; Borrielo & 
Salucci 2001; Salucci 2001) and in elliptical galaxies (Kee- 
ton 2001), although the evidence is less direct than in the 
case of dwarf and LSB galaxies. 

Assuming the non-singular isothermal or pseudo- 
isothermal halo models, several authors have tried to find 
scaling laws for halo cores. In a pioncring paper, Kor- 
mendy (1990) derived the halo component parameters of 
some dSph and spiral galaxies, and concluded that halos 



with larger central velocity dispersions, ct, have larger core 
radii, r c , and smaller central densities, p c (see also Burk- 
ert 1995; Salucci & Burkert (2000)). The correlations they 
found have a large scatter. On the other hand, for nor- 
mal spirals, the halo parameters depend strongly on the 
assumed stellar M/L ratios, and for dSph galaxies the ob- 
servational uncertainties are large. 

More recently, Firmani et al. (2000, 2001a), in an effort 
to infer the core halo scaling laws, have analyzed a sam- 
ple of dwarf and LSB galaxies with high-quality rotation 
curves and cluster of galaxies studied with gravitational 
lensing. They used the non-singular isothermal model and 
took into account the effect of gravitational drag by the 
baryon matter on the halo for the LSB galaxies. Firmani 
et al. found that: (i) p c exhibits a large scatter, while a 
correlation with the maximum circular velocity, v max (or 
ct), is not obvious, and (ii) r c and v max tend to correlate 
with a slope smaller than 1. The different mass mod- 
els of observational samples presented by several authors 
lead to a similar conclusion (Cote et al. 2000 for dlrr's; 
Verheijen 1997, and de Blok et al. 2001b for LSB galax- 
ies). All these authors use different M/L ratio assumptions 
(not too important for dwarfs) and apply either the non- 
singular isothermal or the pseudo-isothermal models. For 
the former case, p c ,-i = 0.54p c and r c ._i = 0.7r c , while 
for the latter p c ,-i = 0.5p c and r Ci _i =r c . We estimate 
Pc-i arid r C) _i from the Firmani et al. compilation and 
the works mentioned above. From Verheijen (1997) and 
de Blok et al. (2001b) we use their fits with the equal 
M/L ratio, and for the Cote et al. (2000) data on nearby 
dwarfs we assume h = 0.65. In the range of the obser- 
vations, v max <~ 40 — 250 kms/s, p c -\ oscillate roughly 
between 10 16 and 10 17 h 2 M Q Mpc -3 with the average at 
s=s 3 x 10 16 h 2 M Q Mpc~ 3 . The radius r C) _i increases with 
Vmax but the scatter is very large; for v max =100 kms -1 , 
r C) _i w 3h _1 kpc on average. Note that for the LSB galaxy 
data from Verheijen (1997) and de Block et al. (2001b) 
the inferred central densities and core radii are upper and 
lower limits, respectively. The formation of a central disk 
contracts adiabatically the inner halo. For a typical LSB 
galaxy, the halo core radius could have shrinked by a factor 
of ~ 1.5 after disk formation (Firmani & Avila-Reese 2000; 
Firmani et al. 2001a), increasing the central density by a 
factor of ~ 2. For higher surface brightness galaxies these 
factors are larger. Observations show that the fitted halo 
central density indeed increases with disk surface bright- 
ness (Verheijen 1997; Swaters 1999). Taking into account 
a correction for the further halo contraction, the scatter in 
the plot p c ,-i-v max diminishes and the correlation between 
r Cl _i and v max becomes tighter. 

At the scale of cluster of galaxies, the highest resolution 
mass distribution analyses were done on the basis of ei- 
ther strong gravitational lensing studies or high-resolution 
X-ray studies of the intracluster medium. For the clus- 
ter CL0024+1654 studied with strong gravitational lens- 
ing techniques (Tyson, Kochanski, & deh" Antonio 1998), 
Firmani et al. (2001a) found a central density similar to 
those of their dwarf and LSB galaxy sample. From Fir- 
mani et al., we calculate p c ,-i = 2.6 x 10 16 h 2 M Q Mpc~ 3 
and r C) _i = 45.5h~ 1 kpc (v max =2400kms~ 1 ). Most of the 
measured cluster mass distributions from high-resolution 
X-ray studies of cluster of galaxies (mainly with the Chan- 
dra satellite) are well fitted by both the NFW profile 
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Fig. 10. — Cumulative v max — function of MW- (left panel) and cluster-sized halos (right panel). All Gl (Cll) subhalos with 
v m ax greater than 20 kms -1 (126 kms -1 ) and within a sphere of radius 200 h _1 kpc (1.5h -1 Mpc) were counted. The averaged 
v max — function from satellites of MW and Andromeda is represented by stars (taken from KKVP99) and it is shown in left panel. 
Error bars are just Poisson errors. The number of subhalos in SI halos with a reasonable odm (i.e., with a odm that also satisfies 
other observational constrains) is comparable to the one seen in the CDM prediction. Notice that there is only one subhalo within 
halo C/I3.0 and in the last two v max bins of halo CZlo.o- The presence of few more massive subhalos in SI cluster-sized halos could 
be explained by a longer survival time of these halos with respect to the survival time of subhalos in the CDM halo. 



and the pseudo-isothermal model. From the fitting pa- 
rameters reported for the clusters A1835 (Schmidt, Allen, 
h Fabian 2001), A2390 (Allen, Ettori & Fabian 2001), 
RXJ1347.5-1145 (Allen, Schmidt, & Fabian (2001)), and 
EMSS 1358+62 45 (Arabadjis, Bautz & Garmirc 2002; this 
cluster is at z = 0.328), we calculate p c _i = 2.8 x 10 16 , 
2.0 x 10 16 , 9.1 x 10 16 and > 2.4 x 10 17 h 2 M Mpc~ 3 , 
r Ci _i = 32.5, 37.5, 22.5 and < 6 h kpc, and v max = 
1843, 1785, 2290 and 1008 kms^ 1 , respectively. Note, 
however, that the inclusion of the adiabatic cluster halo 
contraction due to a cD or another central galaxie would 
favor more extended shallow cores. 

The dotted lines in Figure 4 encompass roughly the 2a 
dispersion of the observational inferences of p c ,-i and r c ,-i 
versus v max for the dwarf and LSB galaxies and the cluster 
of galaxies described above. The dashed lines arc lines rep- 
resentative of the average behavior of the data. In Figure 
4 were also plotted the model predictions for different val- 
ues of a dm- At cluster scales the ACDM halos seem to be 
in agreement with observations. However, at galaxy scales 
they are roughly a factor of ten denser at the radius r Cj _i. 
Extrapolating to scales of dwarf galaxies this factor would 
be even larger. On the other hand, a SIDM model with an 
effective gdm as large as 3.0 cm 2 gr~ 1 produces halos with 
a large variety of inner structures. As discussed in §3, 
for these values of <jdm the core collapse is triggered in a 
time scale smaller than the Hubble time. However, we also 
have seen that the core collapse could be delayed by late 
major mergers or even stopped by a continuos injection 



of thermal energy from the hotter particles of the ambi- 
ent in which the halo is embedded. Thus, in a cosmogony 
with a relatively high effective cross section, some halos 
could show shallow inner density profiles (e.g., elliptical 
and dwarf galaxy halos), while those which are isolated 
would be in the core collapse phase showing cuspy inner 
density profiles (e.g., LSB galaxy halos). Nevertheless, the 
major conflict in this case is that substructure is strongly 
erased; for example, for a dm — 3cm 2 gr~, most subha- 
los in cluster- and galaxy-sized halos are evaporated, in 
complete disagreement with observations (Fig. 10). 

The SIDM simulations with (a ,a) = (0.1,0.0) and 
(0.5, 0.0) show that at galaxy scales the central densities 
of the halos agree roughly with observations, the latter 
case producing probably slightly shallower halos than ob- 
served. However, at the cluster scales the halos in both 
cases are too shallow w.r.t. observations. We are then left 
with SIDM models in which udm varies as the inverse of 
V12. For the cases explored here, (<Jo, a ) — (0.5,1.0) and 
(1.0,1.0), halos at galaxy and cluster scales fall within the 
obsevational range, and p c and the halo soft core fraction 
are approximately constant with scale (v mlix ). Unfortu- 
nately, we can not explore the inner density profiles of the 
(sub) halos at the scales of dwarf galaxies because of they 
are resolved with too few particles. Extrapolating, it seems 
that halos at the scales of dwarf galaxies would remain as 
shallow as the larger ones. A weak hint that this could be 
happening comes from meassuring the cm concentration 
parameter for subhalos of MW-sized halos. This parame- 
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ter characterizes the shape of the density profiles at radii 
<~ 0.03 — 0.1K vir . In the central panels of Figure 7 one 
sees that cm for subhalos of halos Gl .5/v, G2 . 5 / v , and 
Gli.o/y is on average lower than that corresponding to 
CDM subhalos. To a first approximation, we expect p c to 
be independent of v max when a dm oc l/vi2 because N co i 
(eq. 4) or the collisional probability P do not depend in 
this case on vi 2 . However, N co i or P also depend on the 
local density and this is typically larger for smaller halos. 

The satellite cumulative v max — function of the halo Gl 
with ((T , a) = (0.5,1.0) is slightly higher than the one 
of the corresponding CDM halo at the low velocity limit, 
and then drops by a big factor, giving good agreement 
with observations for v max ;> 30kms _1 (Fig. 10). For 
the (do, a) — (1.0, 1.0) case, the function lies above of the 
CDM prediction up to v max =50 kms -1 . As explained in 
§4, this could be due to the softer inner gravitational po- 
tential of the SIDM halos with respect to the CDM ones. 
For the explored range of the SI parameter space, it seems 
that the main influence on the number and structure of 
subhalos comes from the scattering properties of the host 
halo rather than those intrinsic to subhalos (collisions only 
between subhalo particles). In §3.2 we have seen indeed 
how structures of subhalos with similar a dm differ because 
the effective a dm of the host halos are different. 

In conclusion, among the SIDM models studied in this 
paper, those with odm inversely proportional to V12 are 
able to produce halo cores shallower than CDM and with 
central densities roughly constant with v max . However, the 
substructure problem for MW-sized halos is not solved. 
Reionization has been proposed as a mechanism able to 
inhibit the formation of dwarf galaxies within the large 
population of predicted small CDM halos (Bullock et al. 
2000), in such a way that the substructure excess would 
not be longer a critical issue. Nevertheless, a potential 
problem with the Tully-Fisher relation of dwarf galaxies 
would remain (Avila-Rcese et al. 2001). 

Recently, HO carefully analized several observational 
constraints for the SIDM model. Their study is mostly 
based on analytical estimates calibrated on basis of pre- 
vious numerical results. These estimates agree with the 
results obtained here for SIDM models with a constant 
odm and with a a dm inversely proportional to vi 2 . If we 
add the constrain given by the satellite overabundance in 
MW-sized halos, then the region of allowed (<7o, a) values 
becomes smaller than in HO since models with a <J 1 also 
would be ruled out. However, as mentioned above, this 
constrain is weaker than the others because the satellite 
excess problem could be solved invoking reionization. 

HO have introduced one extra constrain, namely the 
formation of supermassive black holes (SMBH) due to SI 
(Ostrikcr 2000). According to their estimates and using 
the observed mass of the central SMBH of our Galaxy and 
the lack of a SMBH in the bulge- less galaxy M33, they 
conclude that <jdm should be smaller than 0.02 cm 2 gr _1 
for any value of a < 1. A critical assumption of this con- 
strain is that the dark halo initially has the typical cuspy 
CDM density profile, i.e., the establishment of the soft 
core by SI occurs only at time scales close to a Hubble 
time. In fact, the mass of the SMBH formed by accretion 
of SIDM is extremely sensitive to the inner halo density 
profile slope. For a slope shallower than —1 this mass is 



much less than 10 3 M Q for a < 1 (see Fig. 1 in HO). Our 
simulations show that the evolving time scales under SI 
for the cases analyzed here are much less than a Hubble 
time. The core expansion phase actually occurs in time 
scales close to the halo dynamical time (Figs. 1 and 5) for 
o dm ~ 1 cm gr . The calculations of halo evolution un- 
der SI of Firmani et al. (2001b) also show that the inner 
halo density profile flattens quickly, as the halo grows and 
virializes. So, accretion of SIDM could not be the relevant 
mechanism for the formation of SMBH and, therefore, the 
strong constraint of HO might not be valid. For the for- 
mation of SMBH other alternative mechanims have been 
proposed (e.g., Granato et al. 2001). 

6. CONCLUSIONS 

A set of high-resolution N-body simulations has been 
run to study the structure and substructure of cluster- and 
MW-sized halos in a SI ACDM cosmology with cross sec- 
tions odm = co(l/vioo) a , where o$ is in units of cm gr , 
and vioo is the relative velocity in units of 100 kms -1 . Our 
main conclusions follows: 

-For low values of a dm, when N co i <; 2 — 5 at z = 0, and 
for a given mass, the lower the a dm, the smaller the core 
radii and the higher the central densities. This behavior 
may be reversed for high values of a dm- In this last case, 
the evolution of halos under SI is fast and the core col- 
lapse phase could begin well before z = 0. We analysed 2 
cluster- and 3 MW-sized halos with (<ro, a) = (3.0, 0.0) and 
different mass assembly histories. We have found that they 
go through the core collapse phase, although this phase 
can be delayed significantly for some halos by both dy- 
namical heating due to major mergers (but not smooth 
mass accretion) and the evaporation of (sub)halo particles 
interacting with the hot particles of the host halo. 

-The SIDM halos expand their cores due to a heat inflow 
from the hotter surroundings. Since the core heat capacity, 
C, is positive, this process leads to the isothermalization of 
the core. After the core maximum expansion, C becomes 
negative down to smaller radii. Then the gravothermal in- 
stability triggers and the core collapses. The system moves 
from a minimum to a maximum entropy state. If the halo 
periphery is heated and the temperature becomes higher 
than in the core, the heat inflows to the core and the core 
expands until the overall halo isothcrmalizcs. This pro- 
cess delays the core collapse, but will hardly reverse the 
gravothermal instability to a runaway core expansion. 

-Present day halos with N co i <J 2 — 5 arc still expand- 
ing their cores or have just had the gravothermal catas- 
trophe triggered. The inner density profiles when isolated 
(host) are flattened by SI only in the innermost regions, 
r < 0.05R„, r . When (Jdm is constant, massive halos are 
suposed to be the most affected by SI (N co i oc pv max ). 
For typical cluster-sized halos to have N co i <; 2 — 5, 
a dm should be smaller than w 0.5 cm 2 gr _1 . For a dm 
0.5 cm 2 gr _1 , the lower the v max , the higher the p c . On the 
other hand, when ct£>m ^ V v i2> this trend is much weaker: 
p c remains almost constant with v max . We compare the 
predicted core densities and radii as well as the scaling laws 
of the halo cores with those inferred from the observations 
of dwarf and LSB galaxies and cluster of galaxies. The 
SIDM models with gdm = 0.5 — 1.0 (1/vioo) cm 2 gr _1 are 
favored. With these values of a dm, N co i is smaller than 
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2-5 from galaxy- to cluster-sized halos. 

-We find that the structure and population of subhalos 
are determined by the interaction between the hot host 
halo particles and the cooler subhalo particles rather than 
by internal processes in the subhalos. Our simulations 
show that overall SIDM subhalos become puffier than their 
CDM counterparts and that they are completely evapo- 
rated in the limit of high a dm- For all the (<r ,a) values 
used here, we do not see signs of a core collapse in sub- 
halos, but we warn that this result can be masked by the 
poor mass resolution of the simulations at the subhalos 
scales. 

-The number of subhalos within MW-sized halos is 
largely suppresed for (<r ,a) = (3.0,0.0). The cumulative 
v max — functions of models (<tq, a) = (0.1,0.0), (0.5,0.0) and 
(0.5,1.0) agree roughly with the observations for v max J> 
30kms _1 , but at the low v max end they lie slightly above 
the corresponding CDM one. For (ao,a) = (1.0,1.0), this 
function is higher than the CDM at all v max . Subhalos 
may survive a longer time in SIDM halos than in CDM 
ones because the inner tidal force in the former is weaker 
than in the latter. Our results show that the range of 
((To, a) values in which the SIDM scenario is interesting 
does not solve the potential problem of the excess of sub- 
structure in MW-sized halos. 

- The SIDM halos have more angular momentum at a 
given mass shell than their CDM counterparts. However, 
this extra angular momentum is small and does not solve 
the conflict between the measured angular momentum dis- 
tribution in CDM halos and that infered from observations 
for dwarf galaxies. On the other hand, the SIDM core ha- 



los are indeed rounder than their CDM counterparts, but 
they are not spherical. This roundness is not only seen in 
the core, it can extend all the way to the virial radius of 
the halo. 

We conclude that the problem of cuspy halos in 
the ACDM cosmology may be solved by introducing a 
modest dark matter cross section inversely proportional 
to the relative velocity of interaction, <jdm ~ 0.5 — 
1.0(l/vioo) cm 2 gr _1 . On the other hand, our simulations 
show that for these values of a dm, the number of subhalos 
in MW-sized halos is close to that of the CDM counter- 
parts. Our simulations also show that SI flattens the inner 
density profiles of growing CDM halos since early epochs 
in such a way that the grow of super massive black holes by 
accretion of SIDM is not an efficient process. A remaining 
interesting question is whether the CDM candidate parti- 
cles may be predicted with the cross sections after their 
decoupling prefered by our results. 
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